Three-dimensional rotation of even-even triaxial nuclei 
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With the self-consistent three-dimensional cranked Hartree-Fock-Bogoliubov (3d-cranked HFB) 
method, various types of rotational motion near the yrast line are investigated in an even-even 
nucleus in the A ~ 130 mass region (sl^Cere). The possibilities of chiral rotations, tilted-rotations, 
and dynamical aspects of these rotations are discussed through the analysis of the 3d-cranked HFB 
solutions. Although a stable planar solution of the chiral rotation is obtained, an aplanar chiral 
configuration is found to be unstable when triaxial deformation is treated self-consistently. 
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Nuclei in the A ~ 130 region are now of special interest 
because it is suggested that a new form of nuclear rota- 
tion can be expected: a chiral rotation 0, |2| ■ This three- 
dimensional rotation was predicted at first in odd-odd 
systems, such as gg^Prys, whose total angular momen- 
tum is built by a particle and a hole in valence orbitals 
in addition to a triaxial rotor, under an assumption of 
the irrotational-flow moment of inertia. Their angular 
momentum vectors point along the shortest (s), longest 
(I) and intermediate (i) axes of the triaxially deformed 
nucleus (7 ~ 30°) jl^, respectively. When the collec- 
tive angular momentum is absent, the total angular mo- 
mentum consists only of the single-particle spins. The 
corresponding state is called "planar" because the total 
angular momentum is in the s — l plane. The planar states 
correspond to "bandheads" of the chiral rotational bands. 
Rotational members built upon the planar state are clas- 
sified as "aplanar" states because these states have the 
genuine chiral configuration. Possible experimental ev- 
idence of the nuclear chirality was reported in iV = 75 
odd-odd isotones in the same mass region through 
findings of pairs of near-degenerate A/ — 1 bands (so- 
called "chiral doublets" ) . 

It is an interesting question whether or not even- 
even nuclei can have the chiral configuration after break- 
ing (Cooper) pairs. To realise the nuclear chirality in 
even-even systems, at least four quasi-particles must be 
excited; two neutrons and two protons. This multi- 
quasiparticle excited state is expected at high spin where 
rotation-alignment occurs due to the Coriolis force. It is 
likely that such a "simultaneous" aligned state would be 
seen at higher spin, compared with the odd-odd chiral 
systems. Experimentally, a candidate for chiral doublets 
in even-even nuclei is suggested in Ip^Ndye H, but there 
is as yet no detailed theoretical analysis about the pos- 
sibility of chiral solutions in even-even nuclei, except a 
brief conference report Q. 

The original model of the nuclear chirality is based 
on the particle (hole)-plus- rotor model which admits the 
presence of a macroscopic rotor with solid triaxial defor- 
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mation (7 ~ 30°). Although there is no concept of a rotor 
in fully microscopic models like the HFB approach, the 
rotor spin can be treated as the collective part of angular 
momentum in the cranking model, which can handle the 
collective and single-particle angular momentum on the 
same footing. An analysis of the chiral configuration in 
odd-odd systems was carried out through the Tilted- Axis 
Cranking model (TAC) 2]. Although each microscopic 
calculation (TAC) is made for given parameters of shape 
and pairing gaps as well as rotational frequency, these 
parameters are carefully determined by a macroscopic + 
microscopic approach (the Strutinski method) in terms 
of the energy minimisation. 

If nuclear chirality is produced in even-even systems by 
four or more quasiparticle excitations, where the quasi- 
particles occupy orbitals with very different orientation, 
it is expected that the collective and individual (quasi- 
)particle motions would be strongly coupled through the 
Coriolis force. Besides, the degree of triaxial deformation 
can be susceptible to the influence coming from rotation- 
alignments at high spin, especially in those nuclei whose 
Fermi level is in the upper shell (i.e., high-fJ Nilsson or- 
bits) ieh. Therefore, if we attempt to describe the 
above physical situations, purely microscopic and self- 
consistent treatments should be made with respect to 
nuclear shape and pairing. In other words, for the aim 
of demonstrating the chirality in even-even systems, it 
is necessary to show the emergence of a genuine three- 
dimensional rotation together with the presence of sub- 
stantial triaxial deformation. Furthermore, compared to 
odd-odd systems, it is numerically advantageous for us to 
study even-even systems because the treatment of pairing 
is more straightforward. 

In this paper, we report a fully microscopic and self- 
consistent study of three-dimensional nuclear rotations 
in even-even nuclei in the A ~ 130 mass region near the 
yrast line [T^ , by means of the three-dimensional cranked 
Hartree-Fock-Bogoliubov method (3d-cranked HFB) 

First of all, let us review the method. The HFB 
equation is self-consistently solved with the pairing-plus- 
Q-Q force under constraints on particle numbers (pro- 
ton and neutron), and angular momentum (length and 
orientation). There are six conditions relevant to the 
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latter constraints. They are, (Ji) = J cos 6', (72) = 
Jsin6'sin(/), (J3) = J sin 6* cos 0, and (Bi) — {i = 
1, 2, 3), with Bi = i (^Qjk + Qkj^ {hi, k: cychc and each 

index runs from 1 to 3). The last three conditions are 
necessary when self-consistent calculations are performed 
in shape, so as to fix the intrinsic coordinate system set 
along the principal axes of the quadrupole tensor at rest 
(i.e., J — 0) "9]. Note that azimuthal (</>) and polar {9) 
angles are measured, respectively, from the 1-3 plane and 
the 1-axis. The nucleus H'^CejQ is chosen for this study 
because of the following three reasons: (i) it is in the 
A ~ 130 region where chiral doublets are seen, (ii) it is 
known to be 7-soft in its ground state ^''^d (iii) it 
has its Fermi levels for protons and neutrons located in 
the high-j intruder orbitals (hii/2) with low- and high- 
components, respectively. According to Baranger and 
Kumar two major shells of the spherical Nilsson ba- 
sis are needed for each sector of isospin (N=4 and 5, for 
proton and neutron), in addition to intruder orbits (ii3/2 
both for proton and neutron). The force strengths for 
the quadrupole and pairing parts are determined through 
the (no-cranking) Nilsson -t- BCS model with the defor- 
mation parameters {(3,j) = (0.162, 0.0°) and gap ener- 
gies (Ap,A„) = (1.069,0.881) [MeV] O. The corre- 
sponding wave function is used as an initial state for the 
self-consistent calculations at J = 0. The HFB equa- 
tion for an arbitrary spin (J) is consecutively solved by 
means of the method of steepest descent through the 
principal- axis cranking procedure {6 = 0°) with a step 
A J = 0.1 ?i until high spin where a four quasi-particle 
state appears with significant triaxial deformation. The 
pairing gaps and quadrupole deformation parameters are 
calculated in a fully self-consistent manner. For in- 
stance, they are calculated to be {(3,j) — (0.132,25°) 
and (Ap,A„) = (0.00,0.18) [MeV] at J = 30;i. The 
self-consistent tilted-axis cranking calculations for J start 
with this principal-axis cranked HFB state as the initial 
state. We confirmed that, similarly to the principal-axis 
cranking, the steepest descent method may be employed 
to solve the HFB equation for given (j) along the direction 
of e from 0° to 90° with a step AO = 0.5°. (j) is given 
every 1° from 0° to 90°. In this way, physical quantities, 
such as the energy surface E'^ {(j>9), quadrupole deforma- 
tions (3''{(t)6), "f^ [(j)9) and pairing gaps A^(06') (r = p, n), 
are calculated for the given spin J as a function of tilt 
angles {9,(1)). 

At first, let us see the energy surface at J = 26?i 
(FigHJ). There is a minimum at (6*, </>) ~ (90°, 45°). 
Deformations are calculated to be (/3, 7) = (0.16,17°). 
Major components of the single-particle spins come from 
both proton and neutron hii/2 orbitals, and they are re- 
spectively calculated as {ji,j2,j3) = (0.0,9.3,2.0) and 
(0.2,2.5,9.6). (The unit for angular momentum is h here 
and hereafter.) The 1-, 2-, and 3-axes correspond to 
the i—, s— , and axes, respectively, in accordance with 
the 7 value. The proton and neutron hii/2 orbitals 
thus point along the s— and /-axes, respectively. Col- 



lective rotation around the z-axis turns out not to be 
present. This situation implies that the corresponding 
states are the "planar" solutions of the chiral configura- 
tion. (There is another minimum at {9,(1)) ~ (47°, 0°), 
but this state is identical to the minimum we have just 
discussed above, from the analysis of the single-particle 
configuration, shape parameters, etc.) 

The "aplanar" chiral configuration would be expected 
at J = 28?i if collective motion could be added around 
the i-axis to the planar states in an "adiabatic" man- 
ner. Adiabaticity here means the same as in the con- 
text of the discussion of rotational bands such as the 
ground state rotational bands (g-bands). Analyses of 
the shape parameters, and of single-particle angular mo- 
menta for hxi/2 protons and neutrons over all the energy 
surface imply that the chiral configuration [l9l | is realised 
for {9,(j)) ~ (70°, 45°) at J = 28?i. The corresponding 
quadrupole deformations are (/3, 7) = (0.13, 32°), and the 
alignments for proton and neutron hj^i/2 orbitals are (1.9, 
9.3, 2.4) and (3.7, 4.3, 11.5), respectively. However, as 
seen in Fig|2 this state is near a saddle point, implying 
that this chiral configuration is unstable. The planar con- 
figuration, that forms an energy minimum at J = 26?i, is 
quite susceptible to the collective rotation, through the 
Coriolis force. 

It is interesting to see what happens when the angular 
momentum is further increased to J = 30?i. Fig|31 shows 
the energy surface at J = 30?i. We have confirmed that 
similar structures to this energy surface are seen up to 
J ~ 387i. There are six shallow minima found in the sur- 
face at J = SOh, including an extremely shallow one. Let 
us label these minima as follows: A: {9,(1)) ~ (58°, 42°), 
Bi: {9,(f>) ~ (28°, 90°), B2: {9,(^) ~ (90°, 25°), Ci: 
9 ~ 0°, C2: {9,^) ^ (90°, 0°), C3: {9,4)) ~ (90°, 90°). 
The groups A, B and C represent a three-dimensional 
rotating solution (which we call "3d-rotation" ) , two- 
dimensional solutions ( "tilted axis rotation" , or "TAR" ) , 
and one-dimensional solutions ( "principal axis rotation" , 
or "PAR" ) , respectively. Let us discuss these minima one 
by one. 

The minima Ci {i = 1, 2, 3): This group is classified as 
one-dimensional rotation, or PAR. The minima are lo- 
cated at {9, (j)) ~ (0°, -), (90°, 0°), and (90°, 90°) for Ci, 
C2, and C3, respectively. The (3 deformation values for 
these minima are almost equal {(3 ~ 0.13). The domi- 
nant components of total angular momentum come from 
proton and neutron hii/2 orbitals in either case. The 7 
deformation values are about 25°, 30°, and —31° for Ci, 
C2, and C3, respectively, so that all of them are strongly 
triaxially deformed. Ci and C2 are identical from the 
analyses of their single-particle configurations, energy, 
deformation, etc. They show a rotation around the i- 
axis while C3 shows a rotation around the s-axis. These 
two types of one-dimensional rotation can be classified 
by analogy with classical mechanics. The moment of in- 
ertia for a rigid body {J^^^) and irrotational flow (J^'" ) 

are then given as J^^ oc 1 — W^/9cos(7 — ^n) and 
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J'^' oc /?2sin2(7 - ^n), respectively (n = 1,2,3). At 
7 — 30°, the i-component of the moment of inertia be- 
comes largest for the irrotational flow, while it is the 
s-component for the rigid body. Ci is thus like "irrota- 
tional flow" and C3 is like "rigid body" rotation. These 
two types are energetically degenerate in the present case, 
so that the irrotational and rigid-body characters coexist 
in this high-spin state. 

The minima Bi {i = 1,2): The attribute of this group 
is two-dimensional rotation, or TAR. Si at {9, </>) ~ 
(28°, 90°) and B2 at (90°, 25°) are likely to be the same 
states (that is, just a rearrangement of the principal 
axes). Their deformations are calculated to be (/3,7) — 
(0.13,29°). The 1-, 2- and 3-axes correspond to s— 
and Z-axes, respectively. Again proton and neutron hii/2 
orbitals play a leading role in building total angular mo- 
mentum. Their components are (6.2,7.0,0.0) for protons 
and (11.2,3.8,0.0) for neutrons. The neutron hii/2 orbits 
have substantial components along the i-axis, making the 
total angular momentum vector for Bi lie in the s — i 
plane. This configuration is different from the "planar" 
chiral configuration having the total angular momentum 
vector in the s — / plane. The energy surface around B 
is shallow along one direction, for example, along (p for 
B2- This structure implies the possibility of fluctuation 
around the TAR. 

The minimum A, {0, (f) ~ (58°, 42°): From the tilt an- 
gles, this "minimum" seems to correspond to the state 
with 3d-rotation, but in fact this is not a 3d-rotating 
state. The quadrupole deformations are calculated to be 
(/3,7) = (0.12,60°), so that the corresponding nuclear 
shape is oblate, that is, axially symmetric rather than 
triaxial. The major contribution to the total angular mo- 
mentum comes from the proton and neutron h.11/2 orbits. 
They are (4.8, 8.3, 1.4) for protons, and (10.1, 4.0, 6.7) for 
neutrons. In this case, the shortest axis (i.e., the symme- 
try axis) corresponds to the 2-axis. The alignment along 
the s-axis comes predominantly from the proton hxi/2 
component, which is consistent with the chiral configu- 
ration in odd-odd systems. However, the neutron hii/2 
component contributes to building angular momentum 
along both the 1- and 3-axes. This is because the princi- 
pal axes of the quadrupole tensor along both of the coor- 
dinates are identical, due to the axial symmetry around 
the s-axis. Therefore, this state is not like the chiral 
rotation, but rather corresponds to TAR of an oblate de- 
formed nucleus. This minimum A is extremely shallow 
(depth is of the order of 10 keV), so that the adjacent 
states are energetically degenerate, implying the pres- 
ence of a "zero mode" . Such a dynamical motion will be 
discussed later. 

The group B contains the lowest minima (hence the 
possible "yrast" state) and the group C has the second 
lowest (E'ox — 10 keV with respect to the energy of Bi). 
At the minimum A, i?ox — 300 keV. As is easily seen 
from the energy surfaces being very shallow, it is more 
natural to think that these states are coupled dynami- 
cally 01 . Namely, the state observed in reality will be 



a superposition of these minima. Nevertheless, it is use- 
ful and meaningful physically to classify these dynamical 
states in terms of several simpler modes. Let us now 
discuss the dynamical aspects of these minima. 

First of all, a fluctuation around Bi should be studied 
because they are the lowest minima. The potential-well, 
for instance, around B2 extends one-dimensionally to- 
ward the direction of (p, indicating a swing motion around 
TAR. This dynamical motion can be called "nutation". 
The nutation is hkely to occur for \(j)\ < 40° with 9 = 90° 
in the case of S2. The signature symmetry, a symme- 
try with respect to 7r-rotation around one of the princi- 
pal axes, say the i-axis, is broken in TAR. However, the 
nutation plays a role in restoring the symmetry dynami- 
cally, as seen in the multi-band crossings in the A ~ 180 
region fl5l|. Another oscillatory motion is expected also 
with respect to 9, although this degree of freedom can 
be treated as a minor correction to the nutation. For 
these dynamical states, the quadrupole deformations are 
calculated to be almost constant and decoupled from the 
tilting degrees of freedom. 

The second type of dynamical motion is expected 
around the minimum A along the shallow plateau. Be- 
cause of the oblate deformation, all the states rotated 
around the s-axis are energetically degenerate. This 
sort of collective rotation is regarded as a zero mode, 
which is closely connected to the restoration of a con- 
tinuous symmetry. In this case, the axial symmetry 
broken by the TAR of an oblate deformed state is re- 
stored by this zero mode. Noting the convention of the 
tilt angles, i.e, (Ji/J, J3/J) = (cos0, sin^cos^), we ro- 
tate the state A around the s-axis (corresponding to the 
2-axis in this case) by a. The rotated state is then 
expressed as 9{a) = arccos( Jf^ cos a/ J — J^sina/J) 
and (j){a) = arccos {esc 6'(a)( Jj^ sina/ J -I- cos a/ J)} . 
Here, ( Jj^, Jj^, J^) denotes the components of total an- 
gular momentum for the minimum A. The trajectory 
{9{a), (j){a)) is drawn in FigO showing that the shallow 
plateau near A corresponds to the zero mode rotation. 
Although the oblate deformed solution collapses beyond 
the plateau, the structure of the energy surface between 
the minima Bi seems to retain the influence from this zero 
mode. If the tunnelling between these minima through 
A has a physical meaning, it should be related to the 
restoration of the broken symmetry, as previously dis- 
cussed by Bonche et al. [T^ . 

These dynamical modes are related to restoration of 
symmetries. In the first case, the nutation of the triaxial- 
deformed state restores the discrete signature symmetry, 
while the zero mode of the oblate-deformed TAR state 
restores the continuous axial symmetry in the second. 
The explicit treatment of these modes should be made 
by means of the generator coordinate method, which it 
is planned to achieve in the future. 

In summary, the self-consistent 3d-cranked HFB equa- 
tion is solved for a 7-soft nucleus ^^''Ce, and several solu- 
tions are obtained for one and two-dimensional rotations 
(PAR, TAR, and planar). All the solutions are built upon 
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four quasi-particle states of hii/2 orbitals for protons and 
neutrons. The chiral configuration is found as a planar 
state (i.e., the "bandhead" of the chiral rotational band), 
but its aplanar partner turns out to be unstable in the 
present case. As the angular momentum increases, the 
corresponding configuration evolves into an oblate de- 
formed state whose rotational motion is TAR-like. This 
oblate state is seen over a relatively wide range of an- 
gular momentum at high spin {30fi ^ J ^ 38fi). As 
a general characteristic, it seems that the rotational mo- 
tion of ^^''Ce at high spin tends to be dynamical and non- 
uniform (nutation). Furthermore, it should be noted that 
the chiral solutions in odd-odd systems are also in quite 
shallow minima in the energy surface (barrier heights 
range from 50 to 100 keV) However, there is an 



argument in the framework of the particle-rotor model 
that the tunnelling between chiral minima is suppressed 
when the mass parameter is large 0- More studies may 
be necessary to understand the dynamical aspects of nu- 
clear chirality, particularly in the context of the coupling 
between deformation and orientation degrees of freedom, 
in theory and hopefully in experiment. 
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FIG. 1: Energy surface at J = 26?i. The contour interval 

in the lower panel is 50 kcV. Discontinuities in the energy 
surface are due to configuration changes, or level crossings. 




FIG. 2: Energy surface at J = 28h. The contour interval in 
the lower panel is 40 keV. 
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FIG. 3: Energy surface at J = 30h. The contour interval in 
the lower panel is 20 keV. The trajectory of the zero mode of 
the oblate state A is schematically displayed by the dashed 
curve without an arrow. The straight line with arrows indi- 
cates the nutation. 



